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Particle acceleration in rotating and shearing jets from AGN 
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Abstract. We model the acceleration of energetic particles due to shear and centrifugal effects in rotating astro- 
physical jets. The appropriate equation describing the diffusive transport of energetic particles in a coUisionless, 
rotating background flow is derived and analytical steady state solutions are discussed. In particular, by consider- 
ing velocity profiles from rigid, over flat to Keplerian rotation, the effects of centrifugal and shear acceleration of 
particles scattered by magnetic inhomogeneities are distinguished. In the case where shear acceleration dominates, 
it is confirmed that power law particle momentum solutions f{p) oc p~^^+°''> exist, if the mean scattering time 
Tc OC p" is an increasing function of momentum. We show that for a more complex interplay between shear and 
centrifugal acceleration, the recovered power law momentum spectra might be significantly steeper but fiatten 
with increasing azimuthal velocity due to the increasing centrifugal effects. The possible relevance of shear and 
centrifugal acceleration for the observed extended emission in AGN is demonstrated for the case of the jet in the 
quasar 3C273. 
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1. Introduction 

Astrophysical jets emerge from a variety of astrophysi- 
cal sources, ranging from protostellar objects, to micro- 
quasars, active galactic nuclei (AGN) and probably 
Gamma Ray Bursts (GRBs). Currently, the observations 
of jets from radio-loud AGN may be counted among the 
most interesting ones as they constitute test laborato- 
ries for studying the spatial structure in relativistic jets. 
Today there is convincing evidence that the central en- 
gine in these AGN is a rotating supermassive black hole 
surrounded by a geometrically thin accretion disk, which 
gives rise to the formation of a pair of relativistic jets. 
The observation of superluminal motions and theoreti- 
cal opacity arguments indicate that the bulk plasma in 
these jets moves at relativistic speeds along the jet axis. 
Models based on the assumption of a one-dimensional ve- 
locity structure have thus allowed useful insights into the 
emission properties of AGN jets. However, as real jets are 
generally expected to show a significant velocity shear per- 
pendicular to their axes, such models may be adequate 
only for a first approximation. In particular, several inde- 
pendent arguments suggest that AGN jet flows might be 
characterized by an additional rotational velocity compo- 
nent: 



(1) The analysis of the jet energetics in extragalactic ra- 
dio sources have revealed a remarkably universal corre- 
lation between a disk luminosity indicator and the bulk 
kinetic power in the jet (Rawlings & Saunders 1991 
Celotti & Fabian 1993) and supported a close link be 



tween jet and disk. The successful application of mod- 
els within the framework of a jet-disk symbiosis (Falcke 
Falcke et al. 



Biermann 1995 



1995) indicates that for 



radio- loud objects the total jet power may approach 1/3 
of the disk luminosity so that a considerable amount of 
accretion energy, and hence rotational energy of the disk 
(cf. virial theorem), is channelled into the jet leading to 
an efficient removal of angular momentum from the disk. 
(2) Observational findings, including helical motion of 
knots or periodic variabilities, seem to provide addi- 
tional eviden ce for intrinsic rotation in AGN jets (cf. 
Biretta |1993| for M87; Ca menzin d & Krockenberger |1992 
for 3C273; Schramm et al. |1993| for 3C345) 
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(3) From a more theoretical point of view, intrinsic jet 
rotation is generally expected in magnetohydrodynamical 
(MHD) models for the forma tion a nd coUimatio n of as tro- 
physical jets (e.g. Begelman 1994; Sauty et al. ^002 ). In 
such models, intrinsic rotation with speeds up to a con- 
siderable fraction of the velocity of light is a natural con- 
sequence of the assumption that the flow is centrifugally 
accelerated from the accretion disk. It should be noted 
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however, that the rotation profile in the jet does not nec- 
essarily have to be disk-like, i.e. the set of available jet 
rotation profiles could be much wider and might include, 
for example, rigid, flat and Keplerian profiles (e.g. Hanasz 
et al. 2000). In particular, rigid rotation inside a well- 
defined light cylinder might be related to foot points of 
the magnetic field lines concentrated near the innermost 
stable orbit (e.g. Camenzind 199(^ ; Fendt 1997a), while 
more generally differential rotation would be intuitively 
expected if there is an intrinsic connection between jet 



motion and the rotating disk (cf. also Fendt 1997b; Lery 
& Frank pOOO| ). 

In this paper, we are interested in the influence of such 
rotation and shear profiles on the acceleration of parti- 
cles in AGN jets. So far, several authors have contributed 
to our understanding of particle acceleration by shear. A 
kinetic analysis was used in the pioneering approach of 
Berezhko (|l98l| ; |l982[ ) and Berezhko & Krymsku ( |l98lD . 
Their results showed that the particle distribution func- 
tion in the steady state might follow a power law if the 
mean interval between two scattering events increases with 
momentum according to a power law. Later on, particle 
acceleration in the diffusion approximation at a gradual 
shear transition in the case of non-relativistic flow veloci- 
ties was studied independently by Earl, Jokipii & MorfiU 
( |1988| ). They re-derived Parker's equation (i.e. the trans- 
port equation including the well-known effects of convec- 
tion, diffusion and adiabatic energy changes), but also 
augmented it with new terms describing the viscous mo- 
mentum transfer and the effects of inertial drifts. Jokipii & 



MorfiU (1990) used a microscopic treatment to analyse the 
non-relativistic particle transport in a moving, scattering 
fluid which undergoes a step-function velocity change in 
the direction normal to the flow. They showed that parti- 
cles may gain energy at a rate proportional to the square 
of the velocity change. Matching conditions in conjunc- 
tion with Monte Carlo simulations for shear discontinuities 



were derived by Jokipu, Kota & MorfiU (|1989| ). The Monte 
Carlo analysis was extended by Ostrowski (1990; 1998), 
who also studied the acceleration at a sharp tangential ve- 
locity discontinuity, including relativistic flow speeds. He 
found that only relativistic flows can provide conditions 
for efficient acceleration, resulting in a very flat particle 
energy spectra which depends only weakly on the scatter- 
ing conditions. The relevance of such a scenario for the 
acceleration of particles at the transition layer between 
AGN jets and their ambient medium was stressed in re- 
cent contri bution s by Ostrowski ( ^0G0| ) and Stawarz & 
Ostrowski ( ^002] ). 

The work on gradual shear acceleration by Earl, Jokipii 



& MorfiU (1988) was successfully extended to the relativis- 
tic regime by Webb (198!;; 1992). Assuming the scattering 
to be strong enough to keep the distribution function al- 
most isotropic in the comoving frame (so that the diffusion 
approximation applies), he derived the general relativistic 
diffusive particle transport equation for both rotating and 
shearing flows. Subsequently, Green's formula for the rela- 
tivistic diffusive particle transport equation was developed 



by Webb, Jokipii & MorfiU ( |1994| ). Applying their results 
to the cosmic ray transport in the galaxy, they found that 
the acceleration of cosmic rays beyond the knee by means 
of galactic rotation might be possible, but not to a suffi- 
cient extent. In its form presented however, it rather re- 
mains an essentially theoretical approach which could not 
be easily related to observations. 



In a previous contribution (Rieger & Mannheim 2001b) 
we developed and applied a new model, that utilizes the 
relativistic transport theory advanced by Webb and that 
permits the analysis of the acceleration of energetic parti- 
cles in rotating and shearing AGN jets. By studying veloc- 
ity profiles known to be typical for such jets, we obtained 
results indicating that the resultant particle energization 
is in general a consequence of both centrifugal and shear 
effects. The formation of power law particle spectra under 
a wide range of conditions reveals the significant potential 
of shear and centrifugal acceleration as a natural explana- 
tion for the origin of the extended, continuous emission re- 
cently obser ved fr om AGN (e.g. Meisenheimer et al. 1997; 
Jester et al. 2001). In the present paper we aim to provide 
a more detailed investigation of this model, including the 
derivation of the relevant particle transport equation (cf. 
appendix B). 

The paper is organised as follows: After a short review 
of the underlying theoretical background in part 2 (see 
also appendix A), we present several applications to rel- 
ativistic jet flows with rigid, Keplerian and flat intrinsic 
rotation profiles in part 3, leaving the detailed derivation 
of the particle distribution function to appendix C. Part 4 
provides a discussion of our results. The observational rel- 
evance of the model presented is pointed out in part 5 with 
reference to recent observations. The paper finally closes 
with a short conclusion. 

2. Model background 

By starting from the relativistic Boltzmann equation and 
by using both the differential moment equations and a 
simple BKG (Bhatnagar, Gross & Krook) time relaxation 
for the scattering term, the relativistic particle transport 
equatio n in t he di ffusion approximation was derived by 
Webb ( |1989| ; |1992| ; cf. also appendix In the underly- 
ing physical picture, it is assumed that scattering of high 
energy particles occurs by small-scale magnetic field ir- 
regularities carried in a coUisionless, systematically mov- 
ing background flow. In each scattering event the parti- 
cle momentum is randomized in direction, but its magni- 
tude p' is assumed to be conserved in the (local) comoving 
flow frame, where the electric field vanishes. Since the rest 
frame of the scattering centres is regarded to be essen- 
tially that of the background flow, particles neither gain 
energy nor momentum merely by virtue of the scattering 
if there is no shear or rotation present and the flow is not 
diverging. However, in the case of a shear in the back- 
ground flow, the particle momentum relative to the flow 
changes for a particle travelling across the shear. As the 
particle momentum in the local flow frame is preserved 
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in the subsequent scattering event, a net increase in par- 
ticle momentum may occur (cf. Jokipii & Morfill |1990D . 
Thus, if rotation and shear is present, high energy parti- 
cles, which do not corotate with the flow, will sample the 
shear flow and may be accelerated by the centrifugal and 
shear effects (cf. Webb et al. |1994[ ). 

For the present application, we consider a rather ide- 
alized (hollow) cylindrical jet model where the plasma 
moves along the z-axis at constant (relativistic) , while 
its velocity component in the plane perpendicular to the 
jet axis is purely azimuthal and characterized by the 
angular frequency Q. Using cylindrical coordinates for 
the position four vector, i.e. x"' = {ct,r,4>, z), the met- 
ric tensor becomes coordinate-dependent (i.e. {gap) = 
diag{—l, 1, r^, 1}) and for the chosen holonomic basis the 
considered flow four velocity may be written in shortened 
notation as 



u" = jfil,0,n/c,vjc), 
Ua = 7/ (-1,0, 17 r^/cjW^/c) 

where the normalization 



7/ = l/v/l-nVVc2-w2/c2 



(1) 

(2) 



(3) 



denotes the Lorentz factor of the flow and where the angu- 
lar frequency may be a function of the radial coordinate, 
i.e. rt = n{r). 



As suggested by Webb et al. (1994), the resulting trans- 
port equation may be cast in a more suitable form if 
one replaces the comoving variable p' by the variable 
$ = ln{H), where H is given by 



H 



p"^ c exp 



dr' 



(4) 



In the case of highly relativistic particles (with p"^ ~ p') 
and using an (isotropic) diffusion coefficient of the form 



(5) 



we finally arrive at the steady state transport equation (cf. 
appendix ^) for the (isotropic) phase space distribution 
function /(r, z,p') 




^^ + (1 

K OZ 



-7/«./0^=--. (6) 



In general, the solution of Eq. (^) can be quite compli- 
cated. However, as we are especially interested in the az- 
imuthal effects of particle acceleration in a rotating flow 
further along the jet axis, it seems sufficient for a first 
approach to search for a z-independent solution of the 
transport equation, i.e. we may be content with an inves- 
tigation of the one-dimensional Green's function, which 



preserves much of the physics involved and which corre- 
sponds to the assumption of a continuous injection along 
the jet (cf. appendix C). In this case the corresponding 
source term becomes 



Q 



— S{t 

Ps 



|5($-$,), 



(7) 



describing mono-energetic injection of particles with mo- 
mentum p' = p'g from a cylindrical surface at r = r^. 
Here, 5 denotes the Dirac delta distribution and the con- 
stant qo is defined as qo — iVs/(8 Tr^p'j^ r^) with the total 
number Ns of injected particles. In order to solve the z- 
independent transport equation, we then apply Fourier 
techniques and consider the Green's solutions satisfying 
homogeneous, i.e. zero Dirichlet boundary conditions. 

3. Applications 

3.1. Rigid rotation - no shear 

In order to gain insight into the particle transport and 
acceleration process, let us first consider the infiuence of 
rigid rotation profiles. One may associate such rotation 
profiles with dynamo action in the inner accretion disk 
around a spinning black hole, which creates a jet magne- 
tosphere filled with disk plasma and rotating with the an- 
gular frequency of its foot points concentrated near the in- 



nermost stable orbit (e.g. Camenzind 1996). The accelera- 
tion of particles due to rigid rotation could represent an ef- 
ficient process for the production of high energy particles. 



This was demonstrated by Rieger & Mannheim (2000) for 
the centrifugal acceleration of charged test particles at the 
base of a rigidly rotating jet magnetosphere for conditions 
assumed to prevail in AGN-type objects. In such cases 
an upper limit for the maximal attainable Lorentz factor 
can be established, which is given by the breakdown of 
the bead-on-the-wire approximation in the vicinity of the 
light cylinder. There are several similarities to the appli- 
cation presented here: If we consider the transport of rela- 
tivistic particles in a rigidly rotating background flow (i.e. 

= Qq = const.) by utilising the above model, shearing is 
absent (i.e. df2/dr = 0) and thus particle energization only 
occurs as a consequence of the interaction with the acceler- 
ating background flow, i.e. due to the presence of centrifu- 
gal effects. It can then be readily shown that for the case 
considered here, Eq. (^) is analogous to the Hamiltonian 
for a bead on a rigidly rotating wire (cf. Webb et al. 1994; 
also Rieger & Mannheim 2001a), i.e. 



H = 



p">c 

If 



(8) 



while the transport equation becomes purely spatial. The 
variable H introduced above, could thus be regarded as 
describing the balance between the centrifugal force and 
the inertia of the particle in the comoving frame. By means 
of Nocther's theorem, H could be shown to be a constant 



of motion, and hence Eq. (C.12) indicates that the ratio 



of particle to flow Lorentz factor is fixed. The comoving 
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particle momentum p' can then be simply expressed as 
a function of the radial coordinate r, i.e. p' = (cf. 

Fig.g). 

We may easily solve the resultant (spatial) transport equa- 
tion analytically, assuming homogeneous boundary condi- 
tions at the jet inner radius and its relevant outer ra- 
dius Tout- In the present case, one may choose the size of 
the inner radius rin to be of the order of the radius of the 
innermost stable orbit around a rotating black hole, while 
the outer radius rout is at any rate physically constrained 
to be smaller than the light cylinder radius (note that for 
relativistic MHD winds, the Alfven point may be close to 
the light cylinder). Results concerning the particle trans- 
port are presented in Figs. |l| and^. They indicate that the 
energetic particles injected at will gain energy by being 
transported in the rigidly rotating background flow, i.e. 
due to the impact of centrifugal acceleration the particle 
momentum increases with position as expected. The pre- 
cise evolution in the immediate vicinity of the light cylin- 
der should however be treated with some caution due to 
the limitations imposed by the applied diffusion approach. 
Fig. |2| reveals a decrease in efficiency if the time be- 
tween two collisions increases with momentum (in case of 
a very weak spatial dependence), i.e. if scattering occurs 
less rapidly for the higher than for the lower energy par- 
ticles. 




0.2 0.4 0.6 0.8 1 

r/ r 



Fig. 2. Spatial distribution N{r)/N{rs) for rigid rotation 
using a different energy dependence of the diffusion coeffi- 
cients, i.e. a = (3 = (solid fine), a = —2, (3 = —0.01 (dot- 
ted line) and a ~ 2, P = —0.01 (dashed line). Boundary 
conditions rin = 0.05 tl, rs = 0.1 tl and rout — 0.999 rL 
has been used for the calculations. 




Fig. 1. Particle momentum p' as a function of the ra- 
dial coordinate in the case of rigid rotation, for parti- 
cles injected at r^ = 0.1 rL with initial Lorentz factor 
js = 10 (solid) and 15 (dashed). Here, rL is defined by 
rL-c(l-z;2/cTVf^o 



3.2. Keplerian rotation - shear dominance 

Let us now consider the acceleration of particles due to 
a rotating background flow with keplerian rotation pro- 



file fl{r) = kr~^^^, where k — \J G M. Again, such flow 
profiles might be associated with jets or disk winds origi- 
nating from the accretion disk around the black hole and 
dragging the Keplerian disk rotation with them. Hence, 
for reasons of self-similarity, keplerian rotation may be in- 
tuitively regarded as one of the most characteristical de- 
scriptions with respect to the velocity profile of intrinsi- 
cally rotating flows. For example, Lery & Frank (2000) 
recently investigated the structure and stability of astro- 
physical jets including Keplerian rotation in the outermost 
part of the outflow and rigid rotation close to its axis (cf. 
also Hanasz et al. 200C ). They also studied the applica- 
tion to non-relativistic outflows from young stellar objects. 
One may thus eventually consider a simple model where 
particles, accelerated in a rigidly rotating flow, are subse- 
quently injected into a Keplerian rotating flow. 
Generally, if we consider Keplerian rotation, both shear 
and centrifugal effects are present. For non-relativistic 
rotation, analytical solutions of the Fourier transformed 
transport equation are given in the appendix C.2. in terms 
of the confluent hypergeomctrical functions. Such solu- 
tions should be appropriate for the outer jet regions. As 
the relative strength of the contribution by shear to that 
of centrifugal energization evolves with r, shear effects will 
eventually dominate over centrifugal acceleration (see dis- 
cussion). This is particularly illustrated in Fig. ^, where 
we have plotted the logarithm of the (normalized) particle 
distribution function / above p' jp'^ > 3 for /3 = (i.e. al- 
lowing no spatial dependence of the diffusion coefficient) 
and different momentum dependence of k. Most interest- 
ingly, well-developed power law momentum spectra are re- 
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Fig. 3. The momentum-dependence of the (normahzed) 
distribution function / for Keplerian rotation using P — 0, 
calculated for a = 0.5, 1.0, 1.5, 2.0, 3.0 at fixed r = 40ri„s- 
Boundary and injection conditions have been specified as 

is 



^in — 10 ?"ms, 

given by as r, 



rout = 1000 Tins, rs 



— 20r,ns, where 



covered, which suggest a spectral slope linearly related to 
the momentum index of the diffusion coefficient, so that 
we have / oc xhis excellently confirms the re - 

sults previously derived by Berezhko & Krymskii ( 198l| ). 
For a coUisionless plasma with a simple shear flow U (y) , 
and by using a simple BGK term, they found the shear 
acceleration to give rise to a power-law momentum spec- 
trum for the steady state comoving particle distribution 
n{r,p') cx p'^ f cx p'~^-^^"\ if the collision time Tc de- 
pends on momentum as Tc oc p'" and a > 0. However, if 
the momentum index a is smaller than zero, i.e. a < 0, 
an exponential spectrum may developed. 

3.3. Flat rotation - interplay between shear and 
centrifugal effects 

In the case of flat rotation (i.e. = ^lor^/r = v^f/r) 
one may easily investigate a more complex interplay be- 
tween shear and centrifugal effects. For in this case, the 
relative strength of shear to centrifugal effects is indepen- 
dent of the radial coordinate and the general solution of 
the fourier-transformed equation could be cast in simple 
analytical terms allowing basic inverse Fourier integration 
to be done, using homogeneous boundary conditions at rin 
and Tout- For the present application, we have considered 
typical jet flows with relativistic Vz/c — 0.95 and different 
azimuthal velocities w^/ (i.e. for a range of 7/ ~ (3 — 4)). 
The corresponding results are plotted in Fig. ^ using a 
constant diffusion coefficient (i.e. a = /3 = 0). Again, the 
calculated distribution functions reveal a powerlaw-type 



momentum dependence, where for low azimuthal veloci- 
ties very steep momentum spectra are recovered, i.e. the 
momentum exponents are found to be near by —8.8 (for 
/30 = 0.075), -6.9 (for = 0.10), -5.1 (for = 0.15) 
and —4.4 (for (3^ = 0.20). The observed flattening of the 
spectra with increasing azimuthal velocities is indicative 
of the increasing impact of centrifugal effects. 




Fig. 4. The momentum-dependence of the (normalized) 
distribution function f{r,p') for flat rotation, calculated 
for nn/rout = 0.02, r^/rout = 0.04 at position r/vout = 
0.2. Chosen azimuthal velocities are Vc/,f /c = 0.075 (short- 
dashed), 0.10 (dotted),0.15 (solid), 0.20 (long-dashed). 



4. Discussion 

The applications presented so far differ in their relative 
contributions from shear and centrifugal effects. We may 
estimate the relative strength of these two more trans- 
parently by comparing their corresponding acceleration 
coefficients: For the shear acceleration coefficient we have 



<P>s 



p/2 Qp 



A 



-ip'^T,T)^iA + a)p'-T 



(9) 



where F deno tes the viscous energization coefficient given 
by Eq. ( |B.10[ ) and cx ro(r)p'" with = X/v'. On the 
other hand, centrifugal acceleration may be described by 
the coefficient 



(10) 



which for highly relativistic particles {p' ~ p"^) and k = 
v'"^ Tc/?> reduces to 



<P>c 



A + a , A 



(11) 
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Hence, for the ratio Cg —< p >c / < p >s of centrifugal to 
shear acceleration we finally arrive at 



2-2 7 nis/^ 



for keplerian rotation 



5f^y/(c^ — vf) for flat rotation 



(12) 



For the case of Keplerian rotation and for the parameters 
chosen above, the relative strength becomes Cs — 0.05, 
i.e. the shear effects essentially dominate over those of 
centrifugal ones. On the other hand, in the case of galactic 
rotation we have Cg ~ (0.3 — 2), which allows a complex 
interplay between centrifugal and shear acceleration with 
increasing azimuthal velocity. 

Our approach utilises the relativistic diffusive parti- 
cle transport theor y as advanced by Webb ( 1989 ) and 
Webb et al. (1994) and thus assumes the diffusion ap- 
proximation to be valid, i.e. the deviation of the particle 
distribution from isotropy to be small. In a strict sense, 
this requires the particle mean free path \cTc\ to be much 
smaller than both, the typical length scale for the evolu- 
tion of the mean (momentum-averaged) distribution func- 
tion and the typical length scale of variation for the back- 
ground flow. Our conclusions are therefore of restricted ap- 
plicability if highly anisotropic distributions are expected 
as, for example, near ultra-relativistic shocks (cf. Kirk & 
Schneider |l^ ; Kirk & Webb |1988D . 

In the application presented here, energy changes as a 
result of radiative (e.g. synchrotron) losses or second-order 
Fermi effects due to stochastic motions of the scattering 
centres have not been considered. One expects the inclu- 
sion of radiative losses to introduce an upper bound to the 
possible particle energy at the point where acceleration is 
balanced by losses, thus leading to a cut-off in the particle 
momentum spectrum. On the other hand, the inclusion of 
second-order Fermi effects would give an additional diffu- 
sion flux in momentum space. It may formally be taken 
into account by a more careful treatment of the scattering 
term in the Boltzmann equation. For non-relativistic jet 
flows, a numerical study of second-order Fermi accelera- 
tion has recently been given by Manolakou et al. ( |1999|) . 
Our present omission of second-order Fermi acceleration 
in the relativistic transport equation appears justifiable 
for the cases where the typical random velocities of the 
scattering centres (as measured in the comoving frame, 
i.e. relative to the flow) are smaller than a product of the 
radial particle mean free path times the rotational flow 
velocity gradient. However, one should note, that estimat- 
ing the effects of second-order Fermi acceleration for the 
case of flat rotation, for example, clearly show them to be 
of increasing relevance for decreasing azimuthal flow ve- 
locities. A more detailed analysis will thus be given in a 
subsequent publication, while the purpose of the present 
model is confined to the analysis of steady-state solutions 
and the essential physical features of shear and centrifugal 
acceleration. 

So far, shear- type acceleration processes in the con- 
text of AGN jets have been investigated by Subramanian 



road suggested by Katz (1991), who considered the parti- 
cle acceleration in a low density corona due to flux tubes 
anchored in a keplerian accretion disk, Subramanian et 
al. (1999) investigated the acceleration of protons driven 
by the shear of the underlying Keplerian accretion disk. 
They demonstrated that the shear acceleration may trans- 
fer the energy required for powering the jet and showed 
the shear to dominate over second-order Fermi accelera- 
tion. However, their model does not deal with the accel- 
eration of particles emitting high energy radiation, but is 
basically confined to the bulk acceleration of the jet flow 
up to asymptotic Lorentz factors of ^ 10. On the other 
hand, Ostrowski (1998, 2000) examined the acceleration 



of cosmic ray particles at a sharp tangential flow discon- 
tinuity. He demonstrated that both, the acceleration to 
very high energies as well as the production of flat parti- 
cle spectra are possible. In order for this model to apply 
efficiently, one requires several conditions to be satisfied, 
including the presence of a relativistic velocity difference 
and a thin (not extended) boundary, as well as a sufficient 
amount of turbulence on both sides of the boundary. It 
seems however that such conditions might be realized, for 
example, at the jet side boundary layer in powerful FR II 
radio sources. 

With respect to shear acceleration, our approach can 
thus be regarded as complementary to the one developed 
by Ostrowski. While our analysis considers the influence of 
a gradual (azimuthal) shear profile in the jet interior, the 
analysis by Ostrowski deals with a sharp tangential shear 
profile at the jet boundary layer. Generally, one expects 
both processes to occur and to contribute to the produc- 
tion of high energy particles, their relative contributions 
being dependent on the specific conditions realized in the 
jet interior and its boundary. 

5. Observational Relevance 

The present results may serve as an instructive example 
revealing the significant potential of shear and centrifugal 
acceleration in the jets of AGN. Since intrinsic jet rota- 
tion is typically expected in the AGN setting, the pro- 
posed mechanism might operate in a natural manner over 
a wide range. This suggests that particle acceleration in a 
rotating and shearing background flow could be of partic- 
ular relevance for an explanation of the continuous optical 
emission observed from several AGN jets (e. g. for 3C273; 
M87; PKS 521-3 65; cf. Meisenheimer et al. |1997| , Jester 
et al. 2001, 2002), as in contrast to shock acceleration, 



this mechanism generally would not be constrained to lo- 
calized regions. Recent observations in fact indicate that 
the radiating particles might be widely distributed so that 
the optical emission from radio jets for example, is not 
confined to bright individual knots (cf. Meisenheimer et 
al. 1996, 1997). In order to explain the apparent absence 



of strong radiative cooling between knots, Meisenheimer, 
Yates & Roser (1997) suggested the operation of an ex- 



et al. (1999) and Ostrowski (1998, 2000): By following the 



tended so-called "jet-like" acceleration mechanism, which 
was thought to be associated with velocity shear and 
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expected to contribute in addition to standard diffusive 
sfiock acceleration. Meanwhile, the need of such an ex- 
tended acceleration mechanism is strengthened by recent, 
high-resolution HST observations of the jet in the 3C273 
(Jester et al. 2001, 2002). Classified as a blazar source, 
3C273 is known as the brightest and nearest (z=0.158) 
quasar, with a jet extending up to several tens of kpc. 
The absence of a detectable counter-jet suggests bulk rel- 
ativistic motion even up to kpc-scales, a conclusion also 
supported by recent models for the X-ray emission ob- 



served with Chandra (cf. Sambruna et al. 2001). The ob- 
servations performed by Jester and collaborators have re- 
vealed an extremely well-collimated jet, where the optical 
spectral index varies only smoothly along the jet indicat- 
ing only smooth variations of physical conditions. Their 
phenomenological analysis yields an optical spectral flux 
index of ^ —(0.65 — 0.8) and provides strong evidence for 
an underlying, universal powerlaw electron distribution, 
which is maintained almost entirely throughout the jet. If 
in the context of shear and centrifugal acceleration, the 
particle energy is locally dissipated by synchrotron radi- 
ation, the spectral emissivity ji, for a power-law particle 
number density distribution n cx p'^ f oc p'^^ is given by 
> oc I/-" with s = (5 - l)/2. Thus, for s = (0.65 - 0.8) 
we require S — (2.3 — 2.6), which for example, may be 
realized for the case of a constant diffusion coefficient and 
high flat rotation or for the case of a simple shear flow 
with moderate momentum-dependent diffusion. 

6. Conclusion 

Observational and theoretical arguments suggest that as- 
trophysical jets should exhibit intrinsic rotation of mate- 
rial perpendicular to the jet axis. Motivated by such ar- 
guments, we have proposed and analysed a basic model 
for the acceleration of energetic particles by centrifugal 
and viscous shear effects which is applicable to relativis- 
tic, intrinsically rotating jet flows. The results obtained 
indicate that particle acceleration in rotating jets might 
represent an active mechanism for the production of the 
synchrotron-emitting high energy particles by giving rise 
to power law particle momentum distributions over a wide 
range of conditions. Due to its anticipated non-localized 
operation, such a mechanism is expected to be of particu- 
lar relevance for the explanation of the recently observed 
extended emission in the jets from AGN. Moreover, as 
a general mechanism this process is able to provide high 
energy particles, which are required as seed particles for 
the (localised) first-order Fermi acceleration widely be- 
lieved to occur in relativistic jets (e.g. Drury 198S ; Kirk 
& Duffy [19991 ). Particularly, in the be ginning the injec- 
tion of seed particles could naturally occur at the base 
of the jet due to centrifugal acceleration (cf. Rieger & 
Mannheim ^000]). 



Acknowledgements. We would like to thank C. Hettlage and 
D. Rieger for a careful reading of the manuscript. Useful dis- 
cussions with S. Jester and M. Ostrowski and helpful comments 



by an anonymous referee are gratefully acknowledged. F.M.R. 
acknowledges support from the German Deutsche Forschungs- 
gemeinschaft (DFG), project number MA 1545/2-2. 



References 

Abramowitz, M., & Stegun, LA. (eds.) 1965, Handbook of 

Mathematical Functions, Dover Publications, New York 
Begelman, M. C. 1994, in: The Nature of Compact Objects 

in Active Galactic Nuclei, eds. A. Robinson, R. Terlevich, 

Cambridge Univ. Press, Cambridge, 361 
Berezhko, E. G. 1981, JETP Letters 33, 399 
Berezhko, E. G. 1982, Sov. Astr. Lett. 8, 403 
Berezhko, E. G., & Krymskii, G. F. 1981, Sov. Astr. Lett. 7, 

352 

Biretta, J. A. 1993, in: Astrophysical Jets, eds. D. Burgarella, 
M. Livio, CP. O'Dea, Cambridge Univ. Press, Cambridge, 
263 

Buchholz, H. 1953, Die konfluente hypergeometrische Funktion 

(Ergebnisse der Mathematik; 2), Springer: Berlin 
Camenzind, M. 1996, in: Solar and Astrophysical 
Magnetohydrodynamic Flows, ed. K.C. Tsinganos, 
Kluwer, Dordrecht, 699 
Camenzind, M., & Krockenberger, M. 1992, A&A 255, 59 
Celotti A., & Fabian A.C. 1993, MNRAS 264, 228 
Drury, L.O'C. 1983, Report on Prog. Phys. 46, 973 
Earl, J. A., Jokipii, J. R., & Morfill, G. 1988, ApJ 331, L91 
Falcke H., & Biermann P.L. 1995, A&A 293, 665 
Falcke H., Malkan M.A., & Biermann P.L. 1995, A&A 298, 375 
Fendt, C. 1997a, A&A 319, 1025 
Fendt, C. 1997b, A&A 323, 999 

Hanasz, M., Sol, H., & Sauty, C. 2000, MNRAS 316, 494 
Jester, S., Roser, H.-J., Meisenheimer, K., et al., 2001, A&A 
373, 447 

Jester, S., Roser, H.-J., Meisenheimer, K., Perley, R. 2002, 

A&A, 385, L27 
Jokipii, J.R., Kota, J., & Morfill, G. 1989, ApJ 345, L67 
Jokipii, J.R., & Morfill, G.E. 1990, ApJ 356, 255 
Katz, J.L 1991, ApJ 367, 407 
Kirk, J.G., & Schneider P. 1987, ApJ 315, 425 
Kirk, J.G., Schlickeiser, R., & Schneider, P. 1988, ApJ 328, 269 
Kirk, J.G., & Webb, G.M. 1988, ApJ 331, 336 
Kirk, J.G., & Duffy, P. 1999, J. Phys. G (Nucl. Part. Phys.) 

25, R163 

Lery, T., & Frank, A. 2000, ApJ 533, 897 

Manolakou, K., Anastasiadis, A., & Vlahos, L., 1999, A&A 
345, 653 

Marcowith, A., Henri, G., Pelletier, G. 1995, MNRAS 277, 681 
Meisenheimer, K., Roser, H.-J., & Schlotelburg, M. 1996, A&A 
307, 61 

Meisenheimer, K., Yates, M. G., & Roser, H.-J. 1997, A&A 
325, 57 

Morse, P.M., & Feshbach, H. 1953, Methods of Theoretical 

Physics, Vol. I, McGraw-Hill Book Company: New York. 
Ostrowski, M. 1990, A&A 238, 435 
Ostrowski, M. 1998, A&A 335, 134 
Ostrowski, M. 2000, MNRAS 312, 579 
Rawlings S., & Saunders R. 1991, Nature 349, 138 
Rieger, F. M., & Mannheim, K. 2000, A&A 353, 473 
Rieger, F. M., & Mannheim, K. 2001a, in: High Energy 
Gamma-Ray Astronomy, eds. F.A. Aharonian & H. Volk, 
AIP ConL Proc. 558, 827 



8 



F. M. Rieger and K. Mannheim: Particle acceleration in rotating and shearing jets from AGN 



Rieger, F. M., & Mannheim, K. 2001b, Proc. of 27th ICRC 

(Hamburg), vol. 7, 2701 
Riffert, H. 1986, ApJ 310, 729 

Sambruna, R.M., Urry, CM., Tavecchio, F., et al. 2001, ApJ 
549, L161 

Sauty C, Tsinganos K., & Trussoni E. 2002, in: Relativistic 

Flows in Astrophysics, ed. A.W. Guthmann et al. 

(Springer: Lecture Notes), in press 
Schramm, K.-J., Borgeest, U., Camenzind, M., et al. 1993, 

A&A 278, 391 
Stawarz., L., & Ostrowski, M. 2002, PASA 19, 22 
Subramanian, P., Becker, P.A., & Kazanas, D. 1999, ApJ 523, 

203 

Webb, G.M. 1985, ApJ 296, 319 
Webb, G. M. 1989, ApJ 340, 1112 

Webb, G. M. 1992, in: Particle acceleration in cosmic plasmas, 

AIP Conf. Proc. 264, New York, 287 
Webb, G. M., Jokipii, J. R., & Morfill, G. E. 1994, ApJ 424, 

158 [WJM 94] 

Wolfram, S. 1996, The Mathematica book. Cambridge Univ. 
Press: Cambridge 



Appendix A: The general steady state transport 
in the diffusion approximation 

In the case where we are concerned with the scattering 
of particles in relativistic bulk flows, it has been typically 
found useful to evaluate the scattering operator in the lo- 
cal Lorentz frame in which the fluid is at rest, i.e. in the so- 



called comoving frame K' (e.g. Webb 1985; Riffert 
Kirk et al. 



198. 



1986 



For in this frame K' , a simple form of 
the scattering operator could be applied if one assumes, 
as in the present approach, the rest frame of the scatter- 
ing centres to be essentially that of the background flow. 
Quantities, which are operated upon the scattering opera- 
tor, e.g. the momentum, might then conveniently be eval- 
uated in this comoving frame, while the time and space 
coordinates are still measured in the laboratory frame K 
characterized by its metric tensor gap- We would like to 
note however, that K' will in general be a non-inertial 
coordinate system (i.e. an accelerated frame) and there- 
fore the related connection coefficients will not vanish. 
With reference to the considered particle transport, the 
covariant form of the Boltzmann equation is thus required, 
which may be achieved by replacing the ordinary (partial) 
space-time derivatives by their covariant derivatives. 
Now, by starting from the relativistic Boltzmann equa- 
tion and by using a perturbation solution of the moment 
equations in the diffusion approximation, i.e. by assuming 
the deviation of the particle distribution from isotropy in 
the comoving frame to be small, Webb ( |1989| ) has derived 
the general equation describing steady state particle trans- 
port in relativistic rotating and shearing flows. Following 
Eq.(4.4) of Webb (1989), the special relativistic diffusive 
particle transport equation for the isotropic, mean scat- 
tering frame distribution (averaged over all momentum 



^ Note, that in the following, quantities which are measured 
relative to K' are labelled with a prime superscript. 



directions) fo{x°',p') =< f > may be written as 
J__d_ 

+V„ (cu"/^ + g") = 



(A.l) 

with a;" the position four vector in the laboratory frame 
K, where the background plasma is in motion with four 
velocity u " , and p' = ml v' the (magnitude of the) parti- 
cle momentum as measured in the local (comoving) fluid 
frame K'. Note, that in order to simplify matters, the sub- 
script is omitted in the following characterisations of the 
isotropic part of the particle distribution function, i.e. / 
now denotes the isotropic part! The total particle energy 
and momentum in the frame K' may be written as 



and 



(A.2) 



respectively, with mo the rest mass of the particle and c 

the speed of light. 

The terms in the flrst line of Eq. (AJ) represent particle 



energy changes due to adiabatic expansion or compression 
of the flow (i.e. the term proportional to the fluid four di- 
vergence Wpu^), due to shear energization (i.e. the term 
involving F) and due to the fact that K' i s an a ccelerated 
frame (i.e. the term cx Uq, cf. also Webb 1985 ). The sec- 
ond line gives the effects of diffusion and convection. In 
Eq. (AJ^), Va denotes the covariant derivative while q" 
denotes the heat flux. This heat flux contains a diffusive 
particle current plus a relativistic heat inertial term oc up 
and is given by 

p' dp' ^ 



q 



(A.3) 



As shown by Webb (1989), Eq. (A.l) could be regarded as 
the relativistic generalization of the non-relativistic par- 
ticle transport equation first derived by Earl, Jokipii & 



MorfiU ( |1988| ). 

The acceleration four vector iir. 



of the comoving (or scat- 



tering) frame in Eqs. (A.l) and (A.3), is defined by 

'Wpu^. (A.4) 



The fluid energization coefficient F in Eq. (A.l) represents 
energy changes due to viscosity. Since the acceleration of 
particles draws energy from the fluid flow field, one ex- 
pects on the other hand the flow to be influenced by the 
presence of these particles. As was shown, for example, 



by Earl, Jokipii & Morfill ( |1988| ) and Katz ( |l99lD , the re- 
sultant dynamical effect on the flow could be modelled by 
means of an (induced) viscosity coefficient. If one considers 
the strong scattering limit, i.e. the case where wTc <C 1, 
with uj the gyrofrequency of the particle in the scatter- 
ing frame (i.e. uj = q B' c/p"^), Tc — 1/vc the mean time 
interval between two scattering events and Vc the colli- 
sion frequency, this fiuid energization coefficient could be 



written as (Webb 1989, Eq. 34) 



F = 



30 



«/3 f 



a/3 



(A.5) 
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where (Tq^ is the (covariant) fluid shear tensor given by 



2 

+ g {9af3 + Up) VsU^ , 



Generally, for a contravariant four vector A" the co- 
variant derivative is given by 



(A.6) 



11/5 dxf3 



(B.6) 



with ga/3 the (covariant) metric tensor. Additionally, for 
the strong scattering limit the spatial diffusion tensor k"^ 
reduces to a simple form given by 



while for the covariant derivative of a covariant four vector 
An one has 



K (g"'^ + W'') , with K^v'^Tc/S, 



(A.7) 



, _ dA^ 



(B.7) 



Hence, for the assumed four velocity Eq. ( p3.3| ), the fluid 
the isotropic diffusion coefficient and v' the comoving four divergence becomes zero, i.e. V/3 — 0, while the 



particle speed. 



Appendix B: Derivation of the steady state 
diffusive particle transport equation in 
cylindrical coordinates 



fluid four acceleration Eq. ( [A.4| ) reduces to 



(B.8 



For the components of the shear tensor Eq. (A.6) we may 
then derive the following relations 



"■no — CTll — <^22 — CaS — Cr02 — cr20 — cr23 — CT32 — 

(7^rVc)|^(l~..Vc^) 



0'12 = cr21 



Consider for the present purpose a cylindrical jet model 
where the plasma moves along the z-axis at constant 
(relativistic) Vz while the perpendicular velocity compo- 
nent is purely azimuthal, i.e. characterized by the angu- 
lar frequency 17, in which case it proves useful to apply uis 
cylindrical coordinates = (ci,r, ^, z). One may then 
choose a set of (non-normalized) holonomic basis vectors The viscous energization coefficient Eq. ( |A.5D then be- 
{ea, a = 0, 1, 2, 3} with = dx/dx", which determine a comes 



ar 



(B.9) 



1— form basis {e"}, known as its dual basis. For cylindrical 
coordinates the metric tensor gap becomes coordinate- 
dependent, i.e. for the covariant metric tensor we have 



15 



dr I 



(B.IO) 



[gap) = diag{-l, l,r'^, 1} , 



(B.l) 



noting that a^^ — ^<^oi, <^^^ = cri2/r'^ and a 



13 _ 



while for the contravariant counterpart is consequently 



As the fluid four divergence vanishes, the first term in 



given by (g"^) = diag{-l, 1, 1}. In particular, aU brackets of Eq. (|A^ becomes zero, while for the second 
partial derivatives of the metric coefficients vanish, except brackets one finds 

for the g22 coefficient, and therefore all connection coeffi- 
cients or Christoffel symbols of second order vanish, except 
for 



9 



22 



T^2 T^2 

^ 21 — J- 12 



(B.2) 



Uaq"^uiq\ (B.ll) 
with the radial particle current = q'' [cf. Eq. ( [A.3| )] 
df 72 r (p'0)2 df ' 



1 _ r 

q =q 



dr 



p' dp' 



(B.12) 



Using the chosen holonomic basis we consider a simple 

plasma flow where the four velocity could be written in (A^) becomes 

coordinate form as u^ea = 7 eo + (7 li/c) 62 -I- (7 Wz/c) 63 
and Uae" = —je° + {-fr'^n/c)e^ + {'-fVz/c)e^, respec- 
tively, or in shortened notation as 



In the steady-state the flfth term in brackets of 



V„ q^ = ~{rqn - « (1 + I'vl/c')^ , 



ld_ 

r dr 



(B.13) 



= j{l,0,n/c,vjc), 
Ua = J (-1,0, ilr'^ /c,vjc) . 

where the normalization 
1 



(B.3) 
(B.4) 



(B.5) 



noting that the third component of the heat flux is given 

by 



/I I 2 2 / 2\ 

^k(1 -f 7 vjc ) — 



(B.14) 



Finally, for the fourth term in brackets of Eq. (A.l) we 
have 



denotes the Lorentz factor of the flow and where the u"Vq/ = (7^2/0) — 



angular frequency may be selected to be a function of the 
radial coordinate, i.e. H, = ^{r). 



dz 



(B.15) 



Now, by collecting together all the relevant terms of 
Eq. (A.l)and by introducing a source term Q (which may 
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depend on r,z and p'), we finally arrive at the relevant 
relativistic steady-state transport equation in cylindrical 
coordinates being appropriate for the considered rotating 
and shearing jet flows: 



J_ d_ 

pl2 Qpl 



P [P ) 3 1 



4 ? 

15 



n 2 / 2x / /4 df 



1 V, - nil + Yvi/c')jr^ + -^{r q'') ^ Q . 
oz oz'^ r or 



(B.16) 



For purely azimuthal, special relativistic flows with 
Vz = 0, the transport equation (B.16) reduces to Eq. (5.2) 
derived in WJM 94. 

As suggested by WJM 94, the derived transport equa- 
tion may be cast in a more suitable form by introducing 
the variable 



$ = ln(i?) 



(B.17) 



replacing the comoving particle momentum variable 
p' . Following WJM 94, we may defin e H such that 
{df/dr)H = —q^/n, with q'^ given by Eq. ( B.12| ), the index 
H denoting a derivative at constant H, i.e. 



H = p'" c cxp - 



dr 



(B.18) 



With respect to a physical interpretation of H we would 
like to note, that in the case of rigid rotation (i.e. fl — 
const.) H could be related to the Hamiltonian for a bead 
on a rig idly rotating wire [cf. WJM 94; see also Eq. (|C.12| )]. 
Hence, by writing /(r, f(r,z,^), the relevant 

derivatives transform like 



^ dfir,z,<i>) 



dr 



df{r, z.p') 



dr 



dr 



dfir,z,p') 



dp' 



q_ 

K 



(B.19) 



and additionally recalling that k = (w'^ ''c)/3 we may ar- 
rive at 



1 dn/dr 
r 



K 



dl 

dr 



J_d_ 

pl2 Qpl 



P (P ) . 



dl 
dr 



r \dr 



+ [3 + a , 
$ \P 



j^nh- (df 



dr 



(B.22) 



where the position and momentum dependence of the col- 
lision time Tc (and thus of the diffusion coefficient) has 
been caught into the definition of the variables a and (3, 
i.e. a and [3 are given by 



a 



d In p' 



and 



91nrc 
91nr 



(B.23) 



respectively. Eq. ( B.22[ ) has been obtained noting that 



c p /{p ) and [dTc/ dp']l Tc — a/p' and using 



1 fdv'' 



dp' 



dp' 



dp' 



(B.24) 



In a similar manner, the terms depending on df /dp' in 
Eq. (B.16) may be rewritten using Eq. ( B.20| ) as 



1 



I p'^Tc 



df{r,z,p') 



dp' 



l^r\ 2,2w.^ .fdnVdfir^^ 



(B.25) 



Likewise, for the terms in Eq. ( |B.16 ) which depend on 
d^f/dp'"^ one finds 



1 / 4 2 
1 I 1 T^^^ V 



' ^ dr J I dp'^ 



5c^ \dr J \ (p'0)2 



a/(r,z,$) 



9$ 



p'^ aV(r,z,$) 



{p 



M\2 



a$2 



(B.26) 



using dp' /dp 



where Eq. (B.21) has been used. 



p'^/p' and noting that p'° c = Now, collecting all relevant terms together the steady 



exp[Jdr' 7^ 51^ r'/c^] exp<i>. As usual, the index <i> in 
Eq. ( [B.19| ) denotes a derivative at constant $. Similarly, 
for the momentum-derivatives we have 



state transport equation (B.16) may finally be rewritten as 



a/(r, z,p') _ a$ df{r, z, $) _ p' df{r, z, $) 



9$ 



{p 



.'0\2 



dp' dp' 
and consequently 

d^f{r,z,p') (p'0)2 -2p'2 a/(r,z,$) 



(B.20) 



dp'- 



(p'0)4 

p'2 d^f{r, z, $) 



(p" 



(B.21) 



Now, using Eq. (B.19) and collecting all terms together 



in the transport equation (B.16) which depend on df /dr 



gj,2 \ r \ p' 



9/ 
dr 

«-2(4 



p 

p'O J a$2 



1^^ df 2 2 / 2',d^f Q 

K dz dz-^ K 



dl 
9$ 



(B.27) 
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with / = /(r, z, $) and where the r-derivatives should be 
understood as derivatives keeping $ constant, i.e. we may 
treat r, z, $ as independent variables. Again, for = 



the partial differential equation (B.27) reduces to the 
steady state transport equation given in WJM 94 [e.g. 
see their Eq. (5.5)]. 

For the present application where we are interested in 
highly relativistic particles with p'" ~ p', the steady-state 
transport equation (B.27) simplifies to 



5V 

Q^2 



1 + 



[3- 



72 n-" 



4 2 

5c2 



IVz df 



dr 



dl 
dr 



Q 



— IJT + (1 + 7 Vz/c ) TT-T = -- 



dz 



dz^ 



(B.28) 



Appendix C: Green's solutions for the steady 
state diffusive particle transport equation 

For general applications one may search for (two- 
dimensional) Green's solutions of the steady state trans- 
port equation (see also WJM 94), i.e. solutions of 
Eq. (B.28) with source term 



Q = qoS{r - rs)6{p' -p's)S{z - z^) , 



(C.l) 



or equivalently (i.e. utilizing the properties of the delta 
function) with source term 



90 



p: 



'f6ir~rs)S{<P-<i>,)S{z~z,), 



(C.2) 



describing monoenergetic injection of particles with mo- 
mentum p' — p'^ at position r — z — Zs- For consis- 



tency, the constant go in Eqs. (C.l) and (C.2) has to be 
defined such, that the relevant expression satisfies the re- 
quirement that 5{r — rs) S{p'—ps ') vanishes unless r — Vg, 
(f) = (f)s, p' = p'g and integrate to unity (or Ng if Ns parti- 
cles are injected) over all space and momentum directions. 
Using cylindrical coordinates this requires 



90 



8 TT^ p'^ 



(C.3) 



Solutions of the transport equation ( B.28| ) may then be 
found by applying Fourier techniques, i.e. by using the 
double Fourier transform defined by 



F{r]fj,,iy)^ dz d$ exp[i (z^ $ + ^ z)] /(r, z, $) , 

(C.4) 

where the inverse Fourier transform is given by 
f{r,z,^)^— dfi / di^ cxp[-i{iy<^ + piz)] 

t:/' J —QO J —00 

xF(r;M,^). (C.5) 



Denoting the Fourier transform of the source term —Q/k 
by Q, we have 

> roc 

dz / (i$ expfi (i^$ + ^z)]-^ 



X 5{r - Ts) (5($ - $s) 6{z - Zs) 



Ks p[ 



- exp[i (i/$s + fJ.Zs)] S{r - r^) , 



(C.6) 



Taki ng th e Fourier transform of the transport equa- 
tion (B.28) one arrives at 



Qj.2 \ r I dr 



5^^ (d^j ([3 + "]- + -) 



IVz . 



K 



F = Q. 



F 



(C.7) 



Let Fg be the (two-dimensional) Green's solution 
of this equation satisfying homogeneous, i.e. zero 
Dirichlet boundary conditions, then Fourier inversion [i.e. 
Eq. ( |C.5| )] yields the required Green's function solution 
fG{r, z,p';rs, Zs,p'g). Generally, the considered Fourier 
techniques yield the Green's function for infinite domains. 
The proper Green's function for bounded domains (i.e. 
for finite z) may be obtained, for example, by the method 
of images (e.g. Morse & Feshbach 1953, pp. 812-816) in 
the case where the boundaries are restricted to straight 
lines in two dimensions or planes in three dimensions. 
The steady state version of the Green's formula (cf. 
WJM 94, Eqs. [3. 24], [7.12]) may then be applied in 
order to arrive at the general solution of the considered 
transport equation. 

Yet, even for the simplification Vz — and 7 = const 
(i.e. using a simple galactic rotation law f2 oc 1/r), the 
Green's functi on sol ution for the steady state transport 
equation Eq. ( B.28 ) is not straightforward to evaluate 
(see WJM 94). However, as we are primarily interested 
in an analysis of the azimuthal effects of particle accel- 
eration in rotating jet fiows, we may be content with 
a z-independent solution of the transport equation, i.e. 
with an investigation of the so-called one-dimensional 
Green's function which preserves much of the physics 
involved (cf. Morse & Feshbach [l95^ , pp. 842-847). This 
(one-dimensional) Green's function may be found by inte- 
grating the two-dimensional ring source Q (which depends 
on the space coordinates r^, z^) from Zj = — cx3 to z^ = 00; 
i.e. the one-dimensional Green's function represents the 
Green's solution for a steady state monoenergetic injec- 
tion of particles from a infinite cylindrical surface parallel 
to the jet axis at radius rg. Hence, by utilizing the fourier 
integral theorem it could be directly shown, that the 
one-dimensional Green's function fc.iD corresponds to a 



z-independent solution of the transport equation (B.28) 
with source term Q{r,p') = go <^(^ ^ i"s) 6{p' — p(,), i.e. 
the (one-dimensional) Green's function we are seeking for 
in the following analysis, represents the solution of the 
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modified transport equation 



1 + 



[3- 



72 n'^r 



4 ? 

5c2 
9o 



(1 



dl 
dr 



= il^<5(r-r,)<5(<i>-$. 



(C.8) 



C.l. Rigid rotation profiles 

In the case of solid body (uniform) rotation shearing in 
the background flow is absent since the fluid moves with- 
out internal distortions. For il = fla ~ const., Eq. 
reduces to the purely spatial transport equation 



5V 
g^2 



1 + P 



(3 + a) 



nlr/c^ 



or 



where the constant fig is defined by 

1^0 



(C.9) 



(C.IO) 



while Qo = -qQS{r- rs)6{^ - $s)/(ksP^), with qo given 
by Eq. (|C.3| ), and where the diffusion coefficient is as- 
sumed to be of the form 







(C.ll) 



with Ko,Po s-iid a,/3 as constants, cf. Eqs. (A. 7), ( B.23| ). 
For rigid rotation Eq. (B.1S|) yields 



Hir,p')^p'^cil-nlryc^-v'Jcy/^ 



7 



(C.12) 



where 7 denotes the Lorentz factor of the flow. As H = 
Hs{rs:P's) is a constant of motion (cf. Noether's theorem, 
see also WJM 94), the particle momentum p' in the co- 
moving frame could be simply expressed as a function of 
the radial coordinate 



p'{r) = mo c t 



m 



ml c4 (1-172 ^2/^2 _^;2/c2) 



(C.13) 



with mo the rest mass of the particle. 

Basically, the solution space of the homogeneous part 



of Eq. (C.9) could be described by a set of two inde- 
pendent solutions, e.g. by the functions j/i(r) and y2{r) 
with Wronskian W{r) = W {yi{r) , y2{r)) — yidy2/dr — 
2/2 dyi/dr, where for an appropriate choice 2/2 (j') = 1- The 
relevant analytical expressions might be directly written 
down for some special cases of interest: a.) For a constant 



diffusion coefficient, i.e. a ~ P — 0, two independent solu- 
tions are given by 2/2 (j') = 1 and 



2/1 (r) - \/l~nlryc^ 



3c2 



In 



cii + ^Ji^nl rye-') 



where for the Wronskian one simply finds 

/ o\ 3/2 



W{r) = -r 



1 - 



(C.14) 



(C.15) 



b.) In the case, where p is negative, the solution yi{r) 
could be expressed in terms of the incomplete Beta func- 
tion (cf. Abramowitz & Stegun 1965| , p. 263), i.e. one fi- 
nally may arrive at the system 2/2 (^) = 1 and 



X B 




for a > —5, and (3 < 
(C.16) 

Note, that now the solutions yi, 1/2 have been defined such 
that the appropriate Wronskian reduces to Eq. (C.15) for 
a = /3 = 0, i.e. we have 

/ AO 9\ (3+a)/2 



Wir) 



1 



(C.17) 



The general (one-dimensional Green's) solution of the in- 
homogeneous differential equation Eq. ( |C.9| ) with monoen- 
ergetic source term Qo defined above could then be written 
as (e.g. Morse & Feshbach |1953|, p. 530) 



f{r,p) = yi{r) 



ki- 



-2/2 (r) 



Qo y2{r) 
W{r) ' 

k2- 



W{r) 



(C.18) 



where fci, are integration constants specified by the 
boundary conditions. 

In the disk-jet scenario the accretion disk is usually as- 
sumed to supply the mass for injection into the jet, thus 
for simplicity one may consider a rather hollow jet struc- 
ture ( cf. M arcowith et al. 1995; Fendt 1997a; Subramanian 
et al. 1999) where the plasma motion in the azimuthal di- 
rection is restricted to a region < r < rout < ''l where 
Tin denotes the jet inner radius, rout the relevant outer 
radius and tl the light cylinder radius. Particles are sup- 
posed to be injected at position rg with initial momentum 
Ps, where ri„ < Tg < rout - By chosing homogeneous bound- 
ary conditions fir — ri„) = and f(r — rout) = 0, the 
integration constants in Eq. ( C.l§| ) are determined by 

fci = -fc2 



and 

2/1 (nn) 

q [2/1 ('"out) - yijrs)] 

1 - 2/l(''out)/2/l(nn) 



(C.19) 
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where 



(C.20) 



with go given by Eq. ( |Cl3[ ), i.e. qo = Ns/ {8tt p'^ Vs). 
Therefore the (one-dimensional) Green's solution may be 
written as 

f{r,p; rs,Ps) yi{r) [ki 6{r - Vs) + ki 9{rs - r) 

-qe{r - Ts)] 

+ [k2 0{r ~ r,) + k2 0{rs - r) + 9yi(r,) e{r - r,)] , 

(C.21) 

where 9{x) denot es th e Heaviside step function. The delta 
function in Eq. (C.2C) and Eq. (C.21) indicates that the 
particle momentum in the comoving frame is directly re- 
lated to the relevant radial position by Eq. (CIS ). In order 
to gain insight into the efficiency of the acceleration pro- 
cess one may introduce a spatial weighting function N(r) 
defined by 

fir,p';r,,p'^)^ N{r)6i<^-<Ps)-= Nir)HsSiH -H,). 

(C.22) 

C.2. Keplerian rotation profiles 

In the case of Keplerian rotating background flow with 
r2(r) = kr~^^'^, k = \fGM generally both, shear and 
centrifugal acceleration, will occur. By applying a simple 



(C.23) 



Fourier transformation (cf. Eqs. |C.4|,| |C75| |), i.e. 

/oo 
d<I> exp[z /, 
-oo 

where the inverse Fourier transform is given by 
1 



/ 



27r 



dt^ exp[— i (f>] F , 



(C.24) 



the transport equation ( |C.8| ) could be written as 



Qj.2 



ai + 02 



OF 
dr 



(C.25) 



where 7(r) is defined by j{r) — 7(r) -^/l — u^/c^, with 
7 the Lorentz factor of the flow, and where Qq denotes 
the Fourier transform of the source term Qq — —qoSij — 
r,)S{<i>-<i>,)/{Ksp',), i.e. 

~ goexp[ii/$s] 



Ks p's 



■ 5{r - Ts) . 



(C.26) 



The abbreviations ai, 02, as, 04 in Eq. (C.25) are defined 

by 

(C.27) 
(C.28) 



ai = (1 + /3), 
0-2 = (3 + a) 
9 



GM 



03 



04 



20 
9 



(c2 

(3 -fa) 



GM 

(c2 - ^2) 



GM 



20 (c2 - t-2) 



(C.29) 
(C.30) 



An analytical evaluation of Eq. (C.25) is rather compli- 
cated. However, a simple set of solutions may be found 
in the case of r being large such that the rotational 
velocity becomes non-relativistic and the approximation 
7(r) = (1 - GM/[{^ - w?)r])-i/2 ^ 1 holds. For, the 
Fourier transformed transport equation then simplifies to 



d^F 

Qj.2 



ai 



a2\dF 1 



dr 



{i V ay, + a^v"^) F ^ Qq 



(C.31) 



which, using the substitution y = 02/^, 02 7^ (i.e. a ^ 
—3), leads to 



y 



d^F 



+ (2 - ai - y) 



OF iv a^ + a^v^ 
dy 



02 



F = 0, (C.32) 



for the homogeneous part of Eq. (C.31). Eq. ( C.32[ ) 
is known in the literature as Kummer's equation (e.g. 
Abramowitz & Stegun 1965| , p. 504). For the general case 
where 02 ^ and (2 — ai) ^ —n, n G Nq, the complete 
solution of this equation, i.e. of the homogeneous part of 
Eq. ( p. 31 ), may be written as 



Pnir, v) = ci /i(r, u) + C2 /2(r, v) 
where the functions /i , /2 are given by 



h{r,v) = M 
f2{r,y) = U 



IV a-j + a^v 02 

,2 - ai, — 

02 r 

i ly a3 + ai 02 

,2 - ai, — 

02 r 



(C.33) 

(C.34) 
(C.35) 



Here, M{a,b,y) and U{a,b,y) denote the confluent hy- 



pergeometric functions (cf. Abramowitz & Stegun 1965, 
pp. 504f; Buchholz |1953| , pp. 1-9), with M{a,b,y) being 
characterized by the series representation 



a (a -1-1) 2 
2!6(6-f 1) ^ 



a{a + l){a + 2) 3 
^3! 5 (6+1) (6 + 2)^ ^ 

while U (a, 6, y) is given by the series 

M{a,b,y) 



(C.36) 



U{a,b,y) 



sinvrft \T{l + a-b)T{b) 

A/(l + a-6,2 



b,y) 



T{a)T{2-b) 



,(C.37) 



with T{x) the Gamma function. 

M{a,b,z) has a simple pole at 6 = —n for a ^ —to or 
for a = —TO if TO > n, and is undefined for b — —n, 
a = —TO and m < n. U{a,b,z), however, is defined even 
for b ±n. 

For the relevant Wronskian one has (e.g. Abramowitz & 
Stegun Il96^, p. 505) 



W{y) = W {M{a,b,y),Uia,b,y)) = - 



my 



r(a) 



(C.38) 
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Thus, for the Wronskian W{r) = W{fi,f2) we find 

-(2-ai) gis/j- 



0.2 



r(a(^.)) 



where 



02 



(C.39) 



(C.40) 



R oo. We may close the path by choosing a rectangular 
contour C which consists of the stated parallel A, the real 
axis, and the outer lines Bi, B2 parallel to the imaginary 
axis at R and —R, respectively. For the relevant ranges of 
a {a —3; 2 — ai 7^ —n) of interest, the integrand has no 
poles within the region bounded by C and thus, by virtue 
of Cauchy's integral theorem, the value of the contour in- 
tegral around C sums to zero. For R ^ 00 the integrals 



For the general solution of the inhomogeneous fouricr "^^^ ^.^^ over B2 vanish as might be verified by using 
equation ( p. 25 ) one finds [cf. Eq. (C.18)] 



F(r, ly) = /i(r, ly) [ki e{r - r,) + fci 0(r, - r) 
+ f2{r,iy) [k2eir-rs) + k2 0{r,-r) 

where q[v) has been defined by 
exp[z V <^s\ 



asymptotic expansion formulas for the integrand. Noting 
that by means of Euler's equation = cosx + i sinx we 
have 



(C.41) 



(C.42) 



f{x'^) cosxdx , 



(C.48) 



and collecting all relevant expressions together, one finally 
may arrive at the integral (2 — ai 7^ — n; a 7^ —3) 

3 



f{r,p';rs,p's) = g{rs,p's,a,(3) exp 
do; cos 



v/20 (3 + a) 



with go given by Eq. ( |C.3D and where fci, k2 are integra- 
tion constants specified by the boundary conditions. For 
homogeneous Dirichlet conditions at the boundaries rjn 
and Tout, i-e. Firin,^) = F(rout, i') = with ri„ < r < Tout 
and Tin < Ts < Tout the integration constants are fixed (for 
each v) and given by 

hif) = -qii^) fi{rin,f) 

/2 (?-out ,>^) hjl's.v) - fl (r-Qut ,v) J2{rs,v) 

h{Tont,t^) /i(nn, v) - /i(?'out, v) h{nn, v) ' 

(C.43) 

/2(nn,I^) 



/i(nn,fs;w^) 



l{Ts,p's,a,l3) exp 



^7v(nn,7'out;W^) 



hN inn, rout; UJ'^] 
for r > rs (C.49) 

3 -f Q , 



do; cos 



«2 • 



(C.44) 



By Fourier inversion, we now may obtain the required 
(one-dimensional) Green's function 



/i(r,c^2)/2(ri„,c.2) 



/lAf(nn,7'out;w2) 



2 /l(rs. Tout; W^) 



= -5- / di^ exp[-ii/$]F(r,z^) , (C.45) 

2 TT 



/iiv(nn,rout;w^) 
for r < Ts 



(C.50) 

where we have introduced the following abbreviations 



with F{r, v) given by Eq. (C.41). As may be obvious from 
the foregoing investigation this fourier inversion is not easy 
to evaluate, not even using numerical methods. However, 
in order to cope with the integration, we may consider the 
following substitution (a 7^ —3) 



(C.46) 



9{rs,p's,a,P) = 



2go?-sv/5(3 + a) ( a2 
Sir Ksp'g r(2 — ai) 



v/20 (3 + a) 



h{ 



(C.51) 
(C.52) 
(C.53) 



M{uj^ +a^,2- ai,a2/r) , 
U{lj^ +a^,2~ai,a2/r), 
/i(nn,w^)/2(rs,cJ^) 

-/2(nn,c.^)/i(rs,c.2), (C.54) 

/2( ''out , 

-/l(w,^')/2(nn,C^'),(C.55) 



for which one finds [cf. Eq. ( |C.40| ) and Eqs. ( |C.27| )-( |C.30| )] 

a(i/) =uj^ + a^, (C.47) 

where a is defined by tr = 3 V3 + a/VSO. with go given by Eq. (C.3). Eq. (C.49) may be evaluated 

Performing the substitution in Eq. (C.45), one arrives at using numerical methods (cf. Wolfram 199f:). /i(r, cj^) and 

an integral with the path of integration A now in the /i(r, o;^) are bounded for uj = while the integrand be- 

complex plane, i.e. parallel to the real axis at a distance haves well enough when cj ^ cx) to allow for a numerical 

3{3 + a)i/[2 ^20 (3 -I- a)], extending from —R to R with evaluation. 
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C.3. Flat rotation profiles 

In the case of flat (galactic- type) rotation with 11 (r) = 
^0 ro I r, the co-operation of centrifugal and shear effects 
could be analysed more directly. Applying Fourier trans- 
formation (cf. Eq. [ C.23 |) the transport equation 
simplifies to 
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d^F a dF b{v) 

1 h — 

where a and h{i') are defined by 
a = + + (3 + 0)72772, 



7 T 



(1 - vl/c') 



(C.56) 

(C.57) 
(C.58) 



with 7 = 1/ ^Jl — rf — vll the Lorentz factor of the flow 
and 7] = V0 / c — il,o rp/c i ts rotational velocity, and where 
Qo is given b y Eq ( |C.26 ). Solutions for the homogeneous 
part of Eq. ( C.56 ) could then be directly written down, 
i.e. one has 



Fi(r,i/) 
Fi(r,i/) 



.i(l-a-V(''-l)'+4bH) 
,i{l-a+^ia-iy+4b{^)) 



with Wronskian W{r, fixed by 
W{r, u) 



Ab{v) 



(C.59) 

(C.60) 
(C.61) 



(C.62) 



The general solution of the inhomogeneous fourier equa- 
tion ( C.5(j ) in the case of homogeneous Dirichlet condi- 
tions at the boundaries rin and Tout, may then be deter- 
mined following the steps given in the previous subsec- 
tion, e.g. see Eqs. (C.41) - ( C.45D . In order to calculate 
the inverse Fourier transform, it proves useful to apply 
the following substitutions 



UJ 



{■i + a)i 




{CM) 



/3 



+ 7^ 



(C.64) 



for which one has 



The complete solution may then be found by proceeding 
as presented in the previous subsection. 



